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Abstract 

A large sample of cosmic ray events collected by the CMS detector is exploited to 
measure the specific energy loss of muons in the lead tungstate (PDWO4) of the elec- 
tromagnetic calorimeter. The measurement spans a momentum range from 5 GeV/ c 
to 1 TeV/ c. The results are consistent with the expectations over the entire range. The 
calorimeter energy scale, set with 120 GeV/ c electrons, is validated down to the sub- 
GeV region using energy deposits, of order 100 MeV, associated with low-momentum 
muons. The muon critical energy in PbW0 4 is measured to be 160+^ ± 8 GeV, m agree- 
ment with expectations. This is the first experimental determination of muon critical 
energy. 



*See Appendix A for the list of collaboration members 
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1 Introduction 



The electromagnetic calorimeter (ECAL) of the Compact Muon Solenoid (CMS) experiment d 
has been designed with a stringent performance goal on energy resolution, driven by the capa- 
bility to detect a Higgs boson decaying into two photons. While waiting to record events from 
beam collisions at the LHC, its outstanding characteristics have been exploited to measure the 
muon stopping power in lead tungstate (PbWC>4) for cosmic ray muon momenta between 5 
and 1000 GeV/c. 

The stopping power for muons in the energy range considered can be conveniently written as 



where E is the total muon energy x is the thickness of the traversed material, commonly mea- 
sured in mass per unit surface, a(E) is the stopping power due to collisions with atomic elec- 
trons, and b(E) is due to radiative processes: bremsstrahlung, direct pair production, and pho- 
tonuclear interactions; a(E) and b(E) are slowly varying functions of E at energies where ra- 
diative contributions are important H2H3J. 

Numerical values of stopping power and related quantities quoted throughout this paper are 
taken from tables in [4]. The tables [see 2 J are obtained from calculations following O includ- 
ing post-Born corrections from [5|, which result in an increase of the critical energy by 6 GeV. 
The theoretical uncertainties are everywhere smaller than the statistical precision of the results 
reported in this paper. The definition of critical energy as the energy at which the average rates 
of energy loss through collision and radiation are equal EUSJ is adopted in this paper. 

The lowest and highest momenta considered in the analysis reported here correspond to re- 
gions dominated by collision and radiative processes, respectively. Comparison of the experi- 
mental values of the stopping power with predictions, extends the validity of the ECAL energy 
scale, set with 120 GeV/ c electron beams, down to energy deposits below 1 GeV. Analysis of 
data in the full momentum range permits a comparison of collision losses with radiative losses, 
thus leading to the first experimental determination of muon critical energy. 

The experimental setup and the characteristics of the data sample are described in Section [2] 
Event selection and measurement of relevant physical quantities are discussed in Section p] 
Possible limitations in the measurement of the muon stopping power are addressed in Sec- 
tion |4] In Section [5j corrections for systematic effects and statistical issues are addressed, and 
results are reported. 



The data set consists of cosmic ray muons traversing the CMS detector located 100 m under- 
ground, collected during the detector commissioning phase. 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam- 
eter. Within the field volume are the silicon pixel and strip tracker, the crystal electromagnetic 
calorimeter and the brass-scintilla tor hadron calorimeter (HCAL). Muons are measured in gas- 
ionization chambers embedded in the steel return yoke. In addition to the barrel and endcap 
detectors, CMS has an extensive forward calorimetry. The silicon pixel and strip detectors pro- 
vide excellent momentum reconstruction, with a 10 % resolution for 1 TeV/ c muons passing 
close to the nominal interaction vertex. 




(1) 



2 Experimental setup and data sample 
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3 Event measurement and selection 



The electromagnetic calorimeter is a hermetic homogeneous scintillator detector made of 
PbW04 crystals. Lead tungstate has a high density (8.3 gem -3 ), a short radiation length, 
(Xo = 0.89 cm) and a small Moliere radius (Rm = 2.0 cm). Crystals are organized in a central 
barrel of 61 200 crystals closed by two end-caps of 7324 crystals each. The light is read out by 
avalanche photodiodes (APD) in the barrel and by vacuum phototriodes (VPT) in the 
end-caps. In the barrel, relevant to the analysis presented in this paper, the individual crystals 
have a truncated-pyramid shape with a mean lateral size in the front section of 2.2 cm and a 
length of 23 cm, which corresponds to 25.8 radiation lengths. The barrel crystals are organized 
in a quasi-projective geometry, where their principal axis is tilted by 3° with respect to the 
nominal interaction vertex, in both the azimuthal and polar angle projections. 

The CMS collaboration conducted a month-long data-taking exercise, known as the Cosmic 
Run At Four Tesla (CRAFT), during October-November 2008 with the goal of commissioning 
the experiment for extended operation (6|. With all installed detector systems participating, 
CMS recorded 270 million cosmic ray triggered events with the solenoid at its nominal axial 
field strength of 3.8 T. An event display of a cosmic ray muon crossing CMS is shown in Fig.[l] 

During data taking at the LHC, the APDs will be operated at a gain of 50. In order to increase 
the sensitivity to low-energy deposits and to reduce the impact of the electronic noise, the APD 
gain was raised to 200 for several of the commissioning runs. In this condition the noise per 
readout channel amounts to about 1 ADC count, which corresponds on average to 9.3 MeV (7|. 

The signal from the APDs, after being amplified and shaped, is digitized at 40 MHz and ten 
pulse samples are recorded for each calorimeter channel. In order to reduce the bandwidth 
required for data acquisition, a reduction procedure is applied. A summary of data reduction 
algorithms is reported in Section |3j for more details the reader is referred to ||8|. 

Because of the angular distribution of the cosmic ray muon flux and the lower sensitivity of the 
ECAL end-caps to the energy released by a minimum ionizing particle, only the barrel part of 
ECAL is used in the present analysis. 

3 Event measurement and selection 

The physical variables to be measured for each muon are the momentum, the path length in 
ECAL, and the energy lost in ECAL. 

The muon momentum is measured from the track fit performed in the inner tracking system (9|. 

The path length in ECAL is estimated by propagating the measured muon track inside the 
calorimeter taking into account bending due to the magnetic field and expected energy loss. 
The average path length in ECAL for the selected muon sample is 22.0 cm, with an rms spread 
of 2.6 cm. 

A dedicated algorithm has been developed to reconstruct the energy deposited by cosmic ray 
muons in ECAL. The main differences with respect to the standard reconstruction algorithm for 
events collected during LHC operation are driven by the requirement of good reconstruction 
and association efficiency down to low energy and for energy deposits associated with muons 
not pointing to the nominal interaction vertex. Online data reduction is based on zero suppres- 
sion (ZS), that is, a readout of channels above the ZS threshold, and on a full readout of selected 
regions (SR) of high interest, as identified by a Selective Readout Processor (SRP), based on the 
information of the Level-1 trigger system of ECAL [8|. For CRAFT runs with APD gain set to 
200, the ZS threshold corresponds to about 20 MeV and the SR threshold to about 170 MeV. 
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Figure 1: An event display of a cosmic ray muon crossing CMS: (a) the full detector and (b) 



detail of the central region. ECAL hits are in magenta, HC AL hits in blue, tracker and muon hits 
in green. The views represent a section of CMS in a plane perpendicular to the beam direction, 
the x axis pointing to the centre of the LHC ring and the y axis pointing vertically upwards. 
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3 Event measurement and selection 



When the condition for SR is met, a matrix of 3 x 3 trigger towers is read out. Each tower is 
comprised of 5 x 5 crystals, the matrix is centered around the tower with energy deposit above 
the SR threshold. The energy reconstruction algorithm identifies clusters starting from a central 
crystal (seed) with an energy deposit of at least 15 ADC counts (139.5 MeV), or from the most 
energetic of a pair of adjacent crystals with at least 5 ADC counts (46.5 MeV) each, the cluster 
energy is determined as the sum of all channels above 2 ADC counts (18.5 MeV) belonging to 
a 5 x 5 matrix of crystals. 

Contiguous clusters are merged and the size of the resulting cluster depends on the impact 
angle of the muon at the calorimeter surface. Clusters in ECAL are then associated with muon 
tracks according to a geometrical match between the track extrapolation to the calorimeter sur- 
face and the centre of gravity of the energy deposit. The energy in the cluster is obtained after 
application of relative calibration constants to individual channels and of an absolute energy 
scale factor. The intercalibration constants are defined to equalize the response of the individ- 
ual channels to the same energy deposits in the corresponding crystals. The absolute energy 
scale was determined from data taken with a 120 GeV/ c electron beam IfTQl . It is defined as the 
response to 120 GeV/ c electrons striking the centre of a 5 x 5 matrix of crystals in a reference 
region of ECAL. With these definitions, the same energy scale is obtained over the entire ECAL 
up to local containment effects. According to Monte Carlo simulation, this containment factor 
corresponds to 97 % of the electron energy for an impact at the point of maximum response of 
the reference crystal matrix. Energy deposits due to non-radiating relativistic particles are fully 
contained in a 5 x 5 crystal matrix, thus the energy scale defined at the test beam is scaled by 
0.97 to measure the ECAL response to muons. Additional corrections, related to containment 
effects due to leakage from the crystals of muon induced secondary particles, are addressed in 
Section |5j 

Only runs in which the APD gain was set to 200 are considered in the present analysis. The 
useful sample of CMS triggers is thus reduced to 8.8 x 10 7 . Furthermore, muons are required 
to cross both the upper and the lower half of ECAL barrel (see Section |5|. The initial sample 
of events considered in the present analysis was thus defined by the combined requirements 
of a muon at the trigger level ITT| , of a single muon track reconstructed by the tracker [9], and 
of an energy deposit associated with the track in both the upper and the lower half of ECAL 
barrel. Moreover, analysis is limited to muons with momentum ranging from 5 to 1000 GeV/ c. 
This interval, which is mainly determined by the size of the event sample, is further divided 
for analysis purposes into 20 equal bins in log(p). For the present analysis, the muon momen- 
tum has to be measured upstream of the energy release in the calorimeter. Thus only energy 
deposits in the bottom half of ECAL are used. The following selections are then applied to the 
initial sample: the uncertainty on muon momentum must be smaller than the bin width; the 
ratio of the energy deposit in the calorimeter to muon momentum must be less than 1 within 
the experimental accuracy; the angle between the muon track at the ECAL surface and the cor- 
responding crystal axis must be smaller than 0.5 radians. These criteria reduce the sample by a 
factor 1.7, the largest reduction coming from the angle selection. Moreover, as will be clarified 
in Section [5j events with energy deposits above 500 MeV in the upper part of the ECAL barrel 
are vetoed. This last selection further reduces the sample by a factor 1.15. The final sample 
consists of 2.5 x 10 5 muons. 



The measured muon spectrum, after the selections, is shown in Fig. 2 'a) The spectrum of 



energy deposited in crystals is shown in Fig. 2 b) 



Figure [3] displays the distributions of AE/ Ax, the measured cluster energy divided by the path 
length in ECAL, for muons with momentum below 10 GeV/ c, where collision losses domi- 
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Figure 2: (a) Momentum spectrum of the muons passing the selections; |(b)| spectrum of the 
reconstructed energy in the lower ECAL hemisphere. A logarithmic binning is used. 



nate |(a)| and for muons with momentum above 300 GeV/ c, where sizable radiation losses are 
expected [(b)] 

4 Measurement of the stopping power 

In making a calorimetric measurement of stopping power, the choice of optimal thickness of the 
sampling material must be the result of a compromise. On the one hand, the thickness should 
be small enough for the variation of primary particle energy to be negligible; on the other, 
the thickness should be large enough to allow a precise measurement of the energy released. 
Furthermore, the correction arising from energy leaking from the measurement volume should 
be kept small. The extent to which the CMS ECAL crystals meet these requirements is discussed 
inlOandO 

4.1 dE/dx approximation 

The stopping power is experimentally determined by the measurement of the energy AE lost 
by a muon when traversing a known thickness Ax, which must be small in order that the 
approximation /(E) = —dE/dx = AE/Ax holds. Quantitatively Af/f, where Af = f(E) — 
f(E — AE), has to be smaller than the desired relative precision on the measurement of /. Since 

then 



(3) 
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4 Measurement of the stopping power 
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Fig ure 3: Measured distributions of AE/ Ax in ECAL; [(a)] for muon momenta below 10 GeV/ c; 
(b) for muon momenta abov e 300 Ge V/ c ; the fraction of events with AE/ Ax > 10 MeV g 

(a) and | (b) | respectively. 
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In the present case, Ax is typically 180 g/cm 2 (about 25 Xo in PbWG^). The derivative of the 
stopping power with respect to E can be written as a sum of the contributions of collision and 
radiation processes, namely 

h i-m)A%)>- 

For the scope of this analysis, b(E) in Eq. 1 can be assumed to be constant, with a value equal 
to that at 1 TeV. Then (d//dE) rad = b = 1.4 x HT 5 cm 2 /g Q. 

This corresponds to a muon radiating 1/400 of the energy when traversing Ax = 180 g/ cm 2 . Its 
average contribution to Af If is then 1 / 400, which is negligible when compared to the precision 
of the measurement. 

The component of the derivative due to collisions is a monotonically decreasing function of E 
in the energy range of the present measurement; its expected value at the lowest momentum 
(5 GeV/ c) is 2.6 x 10~ 5 cm 2 / g [4j. Thus A/// is everywhere small when compared with other 
sources of uncertainty. 

The above arguments apply not only to muons losing an energy equal to the expectation value, 
but also to the case of large fluctuations occurring in radiative processes. In fact /(E) is essen- 
tially a linear function of E in this energy range. 

4.2 Containment effects 

Energy lost by a muon is transferred to energy of secondaries. Since secondaries and their 
daughters may travel a significant distance, the energy lost by a muon traversing a given vol- 
ume is in general different from the energy deposited in that same volume. The energy AE 
considered in the previous section is the energy lost by a muon along a given path, while the 
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first output of the measurement procedure is the energy deposited in a given volume of the 
calorimeter. Taking as a reference volume a cylinder with its axis along the muon direction, 
energy flow across the lateral surface can be made negligible by choosing a sufficiently large 
radius. On the other hand there is always a non-negligible flow of energy carried by secon- 
daries entering and leaving the two ends of the cylinder, since secondaries are produced along 
the entire muon path. In a homogeneous material an equilibrium between the energy flowing 
in through the front surface and the energy flowing out of the rear surface is approximately 
reached when the front surface is preceded by at least an amount of material corresponding 
to the energy-weighted average range of secondaries. In fact, considering a muon entering a 
dense homogeneous material from vacuum, the average energy deposition per unit thickness is 
given by the convolution product of a step function equal to the muon stopping power (muon 
passing from vacuum to homogeneous material) with the function S(t) describing the average 
fractional energy deposition density as a function of depth t for secondaries (i.e. the shower 
depth-profile for high-energy electrons and photons). This convolution product at depth x in 
the material is equal to the muon stopping power times the definite integral from to x of S(f). 
Considering, as an illustration, the non-physical case of a muon whose only interaction with 
matter is photon radiation, at a depth in the dense material at which an electromagnetic shower 
would have deposited 90 % of its energy, the ratio between the average energy deposition den- 
sity and the stopping power would be 0.9. Therefore, considering the spectra of secondaries 
produced by muons, it is found that the equilibrium condition is approached for (5-electrons 
after the muon has traversed a thickness of about 10 g/ cm 2 , while for radiated photons it is 
approached after about 20 radiation lengths. 

In the present case, the situation is well described as a homogeneous region of PbWC>4 preceded 
by order of 10 g/cm 2 of nearby material. Because of the presence of a strong magnetic field, 
material separated by more than a few tens of centimeters from the crystals is less effective in 
producing secondaries flowing into the measurement volume. Thus the crystals diametrically 
opposite to those under study (corresponding to a separation of approximately 2.6 m) have 
little effect on the measured energy loss. Therefore, net containment corrections are expected 
to be small or negligible for collision losses, but sizable for radiative losses. A quantitative 
estimate of these corrections has been obtained by means of a dedicated simulation based on 
the Geant4 package [12], as detailed in Section[5j 



5 Data analysis and experimental results 

This section first addresses corrections to systematic effects; instrumental effects are considered 



m 



5.1 the physics effects identified in SectionW 



5.3 



are considered in 5.2 Statistical analysis of data 



and experimental results are then discussed in 
5.1 Instrumental effects 

The relevant instrumental effects are related to the online data reduction and to the energy 
reconstruction procedure described in Section [3] 

The presence of thresholds introduces a bias in energy reconstruction: noise fluctuations above 
threshold contribute a positive bias, while energy deposits below threshold contribute a nega- 
tive bias. 

When the muon direction is close to the crystal axis, the energy deposited in ECAL is in most 
cases above the SR threshold (90 % of events with the angle between muon and crystal axis 
smaller than 0.1 radians). When this condition is met, the bias arises mainly from effects asso- 
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5 Data analysis and experimental results 



ciated with the clustering threshold and the noise spectrum, signals below threshold making 
a negligible contribution. The probability of noise fluctuation above threshold is measured in 
a dedicated analysis of the noise spectrum. A 14.7 MeV bias is estimated for muons at angles 
smaller than 0.1 radians. At larger muon angles to the crystal axis, the average energy de- 
posit per crystal decreases. Therefore conditions for ZS read-out occur more frequently, with 
the higher ZS threshold reducing the noise bias. Moreover, the probability that the energy de- 
posited in a single crystal is below clustering threshold increases, thus contributing a negative 
bias. Bias is then expected to decrease with increasing angles, but its estimate depends on sev- 
eral contributions, and is less reliable than the estimate at small angles. Both a Monte Carlo 
simulation and the analysis of crystal multiplicity versus angle in experimental data indicate 
a constant positive bias (plateau) at small angles, followed by a monotonic decrease at larger 
angles. 
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Figure 4: Dependence of the raw (dE/dx) on the angle a. between the muon direction and the 
crystal axis, for muon momentum between 5 and 10 GeV/ c. Vertical bars represent statistical 
errors. 



Bias dependence on angle for non-radiating muons has been further studied by analyzing the 
dependence of the raw dE/ dx on angle for muons below 10 GeV/ c in experimental data. The 
result, shown in Fig. |3J is consistent with expectations, although it is not inconsistent with a 
wider plateau, extending up to 0.2 radians, or with a linear dependence of bias on angle with 
no plateau. A bias correction, dependent on angle, is thus applied to the estimated collision 
component of the stopping power, by applying the following procedure. At small angles, the 
bias to be subtracted is 14.7 MeV. At larger angles, the bias to be subtracted is estimated from 
a fit to the data shown in Fig. |3J assuming a plateau followed by a linear decrease with angle. 
The central value of the bias estimate is obtained assuming that the plateau extends up to 
0.1 radians, the hypotheses of no plateau and of a plateau up to 0.2 radians are taken as extreme 
cases in the evaluation of the systematic uncertainty of bias subtraction. The uncertainty in the 
evaluation of the 14.7 MeV bias and the statistical uncertainty of the fit also contribute to the 
systematic uncertainty of the bias correction. The overall systematic uncertainty is estimated 
to be 3.5 MeV, which corresponds to about 1.2% of the average energy deposited by a low- 
momentum muon. 

The bias is estimated to be negligible in the presence of the larger energy deposits associated 



5.2 Containment corrections 
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with radiative processes; it is in any case negligibly small compared to the systematic effects 



discussed in Section 5.2 Therefore, no bias correction is applied to the expected contribution 



to dE/ dx from radiative processes over the entire muon momentum spectrum. 

5.2 Containment corrections 

Energy leakage and the effects of upstream and downstream materials are quantitatively dif- 
ferent for collisions and radiative processes. Therefore, in order to make an experimental de- 
termination of the muon critical energy, losses and compensation for radiative and collision 
processes are treated separately. 

To estimate the difference between energy lost and energy deposited in ECAL, a model based 
on the Geant4 simulation toolkit, with a simplified description of the detector and materials 
surrounding it, was developed. The use of a simplified simulation instead of the full simula- 
tion available in CMS was driven by the high statistics needed for these studies and for those 



reported in Section 5.3 The results of the simulations were also compared with dedicated ana- 
lytical calculations. 

Because of the effect of the intense magnetic field on the trajectories of charged particles, ma- 
terial lying more than a few tens of centimeters away from the crystals under study may be 
neglected. In particular, the upper half of the ECAL has been ignored in the simulation. The 
validity of this approximation has been reinforced by selecting muons that produce a signal 
of less than 500 MeV in the upper part of the ECAL, thus removing most of those emitting a 
hard photon. As a further check, the energy deposited in the lower half of ECAL by muons 
releasing less than 500 MeV in the upper ECAL hemisphere (the selected sample) has been 
compared with the complementary sample and agreement was found within the precision of 
the comparison. In fact, in the muon momentum region between 145 and 1000 GeV/ c, where 
radiative losses are relevant, the weighted average of the ratio of the dE/ dx measured in the 
two samples in seven momentum bins is 0.95 ± 0.09. This result gives no indication of an ef- 
fect of the upper part of ECAL on the energy deposited in the lower part, thus supporting the 
hypothesis that far material can be neglected. However, due to the modest statistical precision 
of this test, it was conservatively decided to activate the veto in the event selection. 

The low-momentum region is of major interest for the present analysis, since the large sam- 
ple size and the small energy loss fluctuations allow precise results to be obtained. According 
to the simulation, the energy carried by secondaries flowing out of the rear detector surface 
is, on average, 3 % of the energy lost in the crystals for a muon with momentum 15 GeV/ c. 
This represents an upper limit to containment corrections in the low-momentum region, where 



collision losses dominate. In fact, as already discussed in Sect. 4.2 the rear leakage can be com- 
pensated by secondaries produced in upstream material entering the front detector surface. A 
comparison of the deposits in the upper and lower hemispheres of the ECAL barrel, for muons 
with momentum between 5 and 10 GeV/c, shows no difference in dE/dx, with a sensitivity 
better than 1 % (7|. Given the similarity of the results, in spite of the asymmetric distribution 
of the outer and inner material around ECAL, the net correction has been neglected, and 1 % 
systematic uncertainty on the null correction for collision losses has been assumed. 

Estimated corrections for energy containment in radiative processes have been derived from 
dedicated simulations of two extreme cases: with no material in front of ECAL and with the 
whole material budget of the tracker concentrated just in front of ECAL. The former case rep- 
resents the upper limit to the containment correction, where only rear losses are considered; 
the latter case gives a lower limit to the correction, as it maximally overestimates the energy 
flow through the front face due to upstream material. In both cases the effects of the upper 
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5 Data analysis and experimental results 



ECAL hemisphere, used as a veto in the analysis, were ignored. The mean value and half of 
the difference between the results of the two simulations are the estimate of the corrections 
for the energy containment in radiative processes and of the associated systematic error. The 
correction for energy leakage in radiative processes is maximal at 1 TeV/ c: at this muon mo- 
mentum it corresponds to (28 ± 5) % of the average energy lost, while at 170 GeV/ c it reduces 
to (14.5 ± 2.5) %, where the error arises from the systematic uncertainty mentioned above. The 
uncertainty is small compared to the statistical uncertainty in the measurement of the stopping 
power in the corresponding region. Given the size of statistical uncertainty and the contribu- 
tion from other sources of systematic errors, the precision of the correction obtained with the 
simplified simulation is sufficient over the entire momentum range. The correction for non- 
containment of radiated secondaries is applied over the whole momentum spectrum to the 
expected contribution of radiative processes to stopping power. 

5.3 Critical energy, energy scale, Cherenkov contribution 

Several possible sources of systematic error have been taken into account in making a compar- 
ison between the energy scale of the ECAL set with a 120 GeV/ c electron beam and the scale 
resulting from the measurement of muon stopping power. 

Four out of 36 supermodules were calibrated with 120 GeV/ c electrons. The energy scale was 
transferred to all supermodules by means of the cosmic ray intercalibration discussed in ||T0|. 
This procedure resulted in a typical single-channel uncertainty of 1.5%, which is dominated 
by statistics, and an overall scale uncertainty whose contribution to the final precision of the 
present measurement is negligible. 

The precise variation in the gain of each APD, when going from a nominal gain of 50 to a 
nominal gain of 200, has been measured by means of the laser monitoring system described 
in (TJ. The resulting uncertainty on the average energy scale is negligible. 

The overall variation of the response to incident electrons with temperature has been measured 
in test beam to be -(3.8 ± 0.4)%°C" 1 . In CMS, the temperature is monitored with a 0.01°C 
precision, and the cooling system has been demonstrated to maintain the temperature stable to 
18 ± 0.1°C JTJIZDr thus uncertainties on temperature corrections are negligible. 

A test of reproducibility of calibration constants was performed by exposing the same module 
to test-beam electrons on two occasions, separated by 30 days. This resulted in no measurable 
offset on the energy scale |(T0j . 

A different pulse height reconstruction procedure is required for cosmic-ray data, which are 
asynchronous with respect to detector clock, compared to collision data, which are 
synchronous, or to test-beam data, which are also asynchronous but were recorded together 
with accurate timing information from trigger. This has been estimated to introduce a 
systematic uncertainty of 1% in the energy scale of cosmic ray data with respect to test beam. 



In Fig. E 'a) the specific muon energy loss resulting after corrections is compared to expectations 



as a function of the muon momentum |@|. Fig. E b) shows the ratio of experimental measure 



ments to expected values. Two regions are indicated in Fig. |5] the expected 68 % probability 
central interval (grey shaded area), and the minimum interval corresponding to 68 % probabil- 
ity (region delimited by continuous curves). The reasons for considering both intervals and the 
procedure adopted for their estimate are summarized below. 

The characteristics of radiative losses, which correspond to rare processes with high energy 
releases, generate probability distribution functions (pdf) of energy loss for single events with 



5.3 Critical energy, energy scale, Cherenkov contribution 
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Figure 5: (a) Muon stopping power measured in PDWO4 (dots) as a function of muon momen- 
tum compared to expectations |4| (continuous black line). The expected contributions from 
collision and radiative processes are plotted as well (red dotted line and blue dashed line re- 



spectively), (b) Ratio of the measured and the expected values of the muon stopping power, 
as a function of muon momentum. As discussed in the text, in both figures, the shaded grey 
area indicates the expected 68 % probability central interval, while the continuous cyan curves 
delimit the minimum interval containing 68 % of the expected results. "Experimental errors" 
are not shown for reasons discussed in the text. 
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5 Data analysis and experimental results 



long tails at high energy, that is: highly positively skewed distributions. In this condition, con- 
vergence to a normal distribution for the average dE/ dx of a sample of N events is slow. The 
skewness of the single event distribution increases with increasing muon momentum, owing 
to the increasing contribution of radiative losses with momentum. The number of events per 
bin in the selected data sample rapidly decreases with momentum for p > 50 GeV/ c (Fig. 2 a)} . 



It is thus expected that, due to the combined effect of increasing skewness of single event dis- 
tribution and decreasing event statistics with increasing momentum, in the higher momentum 
bins the pdf of the sample mean of dE/ dx cannot be approximated by a Gaussian. Moreover, 
when this condition occurs, sample mean and sample variance are highly correlated. It fol- 
lows that the widely used "experimental error", namely the sample rms divided by the square 
root of the sample population, is not a good measure for the consistency of the sample mean 
with the expected value. Therefore the uncertainty of the result has been estimated, consis- 
tently over the entire momentum range covered by the measurement, by means of a numerical 
technique based on the simulation of the expected results with Geant4. Ten thousand identi- 
cal experiments, each with the same statistics per momentum bin as the actual data sample, 
were simulated. The expected pdf of the mean is then obtained for each bin as the distribution 
of the mean values of the stopping power from the different experiments. For each pdf two 
68 % probability intervals for the expected result were derived: the central interval, obtained 
by discarding 16 % of the results on each tail of the pdf, and the interval of minimum width 
containing 68 % of the results; this last interval always contains the most probable value. For 
Gaussian distributions these intervals coincide. In Fig. |5j for the event samples corresponding 
to the bins above 100 GeV/ c, the differences between the two intervals become increasingly 
significant. At lower momenta, where the two intervals tend to coincide, they also correspond 
to the conventional estimate of the experimental error. 

For the high momentum bins, the increasing probability of large fluctuations in radiative en- 
ergy losses, combined with the decreasing size of the event samples, cause the most probable 
values to be systematically lower than the expectation values. This trend is particularly marked 
in the highest momentum bin of Fig. [5j where the expectation value lies outside the 68% prob- 
ability interval. 

The curve 



(dE/dx), 



dE\ „ /dE\ 

dx )coll X \ dx )rad. 



where coll and rad label the predicted energy losses in PbWG*4 due to collisions with atomic 
electrons and radiative processes respectively [4], is fitted to experimental stopping power data 
using a binned maximum likelihood and the pdf described above. The parameters a and /3 
account for the overall normalization of the energy scale and for the relative normalization of 
radiation and collision losses. With the adopted parameterization the /3 parameter, from which 
the critical energy is measured, is independent of the overall energy scale. The fit results in: 



a = 1.004+°;°^ ( sta t.) ± 0.019 (syst.) 
f> = 1.07±g;g| (stat.) ± 0.6 (syst.). 

Adding statistical and systematic contributions in quadrature, it may be concluded from the 
above results that the energy scale is consistent with expectations within a systematic uncer- 
tainty of 1.6 %, with a 1.2 % contribution from the uncertainty on the energy scale dependence 
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on the angle and on the clustering, and a 1.0% contribution from uncertainty in containment 
corrections. This result is mainly determined by the precision of the measurements in the mo- 
mentum region below 20 GeV/ c, where radiation losses are marginal. The typical energy re- 
leased in this region is about 300 MeV. Thus the energy scale set with 120 GeV/ c electrons still 
holds in the sub-GeV region. This, assuming a linear response of the detector to the collected 
light, implies that the relative contribution to collected photons of Cherenkov light and scintil- 
lation are the same for 120 GeV/ c electrons and for muons in the 5-20 GeV/ c region, within the 
1.6 % precision of this measurement (this last statement strictly holds only for muons at small 
angle with respect to the crystal axis). 

From the fit results, a muon critical energy of 160+^ (stat.) ± 8 (syst.) GeV is obtained, m agree- 
ment with the computed value of 169.5 GeV for PbW04 [4]. The systematic uncertainty in- 
cludes a contribution of 4.5 GeV from the uncertainty of the containment corrections, domi- 
nated by the limited knowledge of the correction for radiation losses, and a contribution of 
6 GeV due to the stability of the result against bias subtraction and event selection. For the 
latter, sizable contributions come from the variation of the acceptance of the muon angle with 
respect to the crystal axis, and, to a lesser extent, from the requirement that E/pbe lower than 
1. The contributions from these sources were estimated by varying the muon acceptance angle 
between 0.3 and 0.7 radians, and by removing the requirement that E/pbe lower than 1. 

It might be argued that the use of pdf, calculated a priori in anticipation of a particular outcome, 
would be inappropriate when the fitted parameters differ from the expected values. In order 
to investigate the possible importance of this, fits were performed with alternative pdf sets, 
parametrized as functions of a and /3, with a controlling the pdf energy scale, and /3 controlling 
the relative contributions of radiation and collisions. The changes in the results were found to 
be negligible compared to the quoted error. 

6 Conclusions 

The muon stopping power in PbW04 has been measured over the momentum range from 
5 GeV/ c to 1000 GeV/ c, yielding the first experimental determination of the muon critical 
energy. The value obtained: 160+J (stat.) ± 8.0 (syst.) GeV, agrees with the value calculated for 
this material. 

In the region corresponding to muon momenta less than 20 GeV/ c, where collision losses dom- 
inate, the average energy deposited in the crystals is of order 300 MeV. Thus the agreement 
(to within 1-2 %) between the measured stopping power and the calculated values confirms 
that the energy calibration of the detector, previously determined with 120 GeV/ c electrons, 
remains valid down to the sub-GeV scale. Alternatively, the observed linearity in energy re- 
sponse may be interpreted as demonstrating that the relative contributions of Cherenkov and 
scintillation light to the measured signal are the same for muons of ~10 GeV/ c as for electrons 
of 120 GeV/ c (in the specific case of muons traversing the crystals at small angles with respect 
to the axis, and of electrons striking the crystals at normal incidence). The confirmation of 
the detector energy calibration and linearity of response relies on the validity of the calculated 
stopping power in lead tungstate, in the region dominated by collision losses. On the other 
hand, the value extracted for the critical energy does not depend on the energy calibration, 
although it would be sensitive to a deviation from linearity of the response. 
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6 Conclusions 
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